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Quintessential inflation — in which a single scalar field plays the 
role of the inflaton and quintessence — from a sum of exponentials 



potential V = A (e 5ip + e^A or V = A (e 5v + e v ) or a cosh poten- 



tial V = 2j4cosh(5</?) is considered in the context of five-dimensional 
gravitation with standard model particles confined to our 3-brane. Re- 
heating is accomplished via gravitational particle production and the 
universe undergoes a transition from primordial inflation to radiation 
domination well before big bang nucleosynthesis. The transition to 
an accelerating universe due to quintessence occurs near z ~ 1, as in 



Braneworld quintessential inflation can occur for potentials with or 
without a minimum, and with or without eternal acceleration and an 
event horizon. The low z behavior of the equation of state parameter 
Wff, provides a clear observable signal distinguishing quintessence from 
a cosmological constant. 

1 Introduction 

In five- dimensional gravitation with standard model particles confined to 
our 3-brane, the Friedmann equation is modified PQ [3] at high energies: the 
square H 2 of the Hubble parameter acquires a term quadratic in the energy 
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density, allowing [4] slow-roll inflation to occur for potentials that would be 
too steep to support inflation in the standard Friedmann-Robertson- Walker 
cosmology. In this model, quintessential inflation — in which a single scalar 
field plays the role of the inflaton and quintessence — can occur for a sum 
of exponentials or cosh potential, a type of potential that arises naturally 
in M/string theory for combinations of moduli fields or combinations of the 
dilaton and moduli fields. 

Inflation ends in this model as the quadratic term in energy density in 
H 2 decays to roughly the same order of magnitude as the standard linear 
term (technically, when the slow-roll parameter e = 1). Reheating can be 
accomplished [5] via gravitational particle production at the end of inflation. 
The universe subsequently j2] undergoes a transition from the era dominated 
by the scalar field potential energy to an era of "kination" dominated by 
the scalar field kinetic energy, and then to the standard radiation dominated 
era well before big-bang nucleosynthesis (BBN). With a sum of exponen- 
tials or cosh potential, the universe evolves at this point according to the 
quintessence/cold dark matter (QCDM) model. The transition to an accel- 
erating universe due to quintessence occurs late in the matter dominated era 
near z « 1, as in ACDM. 

We will assume a flat universe after inflation. In the QCDM model, 
the total energy density p = p m + p r + P</> = p c , where p c is the critical 
energy density for a flat universe and p m , p r , and p^ are the energy densities 
in (nonrelativistic) matter, radiation, and the quintessential inflation scalar 
field 0, respectively. Ratios of energy densities to the critical energy density 
will be denoted by Q m = p m j p c , Q r = p r /p c , and fl^ = p^/pc, while ratios 
of present energy densities p m o, p r o, and p^o to the present critical energy 
density p c0 will be denoted by fl m o, fi r o, and fl^o, respectively. 

Using WMAP3 [Bj central values, we will set fi 0o = 0.74, fi r0 = 8.04 x 
1(T 5 , Q m0 = 1 - tttpo - tt r0 « 0.26, and p^J 4 = 2.52 x 10" 3 eV, with the 
present time t = 13.7 Gyr after the big bang. 

For the sum of exponentials potential, we will take a monotonically in- 
creasing function of <p (<p will evolve from <pj 3> 1 at the end of inflation to 
|<Po| ~ 1 today) 



where A and B are positive constants, A > p > 0, <p = <f>/Mp, the (reduced) 
Planck mass M P = 2.44 x 10 18 GeV, and with A = 5, p = ^2 or 1, and 
B = 1 for the simulations. We will also discuss the simple exponential 
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potential V(<p) = Ae Xip considered in Ref. [5] and the cosh potential V(<p) = 
2A cosh(A<^) with A = 5. For all four potentials, V ~ Ae Xip for ip ^> 1 (which 
will be the case during inflation and gravitational particle production). The 
constant A ~ p c o for quintessential inflation (see Table 1). 

Note that the BBN (z ~ 10 9 -10 11 ), cosmic microwave background (CMB) 
(z ~ 10 3 -10 5 ), and large-scale structure (LSS) (z ~ 10-10 4 ) bounds % Z 0.1 
are satisfied in all the simulations below (Figs. [21 El arid [T2"j) and that the 
transition from the era of scalar field dominance to the radiation era occurs 
around z pa 10 15 /A. 

Quintessential inflation in the standard cosmology from a sum of expo- 
nentials potential was considered in Ref. [7] for A = 20 and /i = 0.5 or 
—20. Ref. [S] extended the analysis to the braneworld case where H 2 has 
a quadratic term in energy density, for A ~ 20 and fi ~ —A. Braneworld 
quintessential inflation is also discussed in Ref. [9], for A = 4 and \i = 0.1. 
Only one simulation for the sum of exponentials potential braneworld sce- 
nario is presented, in Ref. [9] for the equation of state parameter w^. We 
find that the case A = 4 and \x = 0.1 violates even the looser CMB bound 
ft>4> < 0.2 (see Fig. [T]), and that the recent average of the scalar field equation 
of state parameter w = —0.68 is too large. 
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Figure 1: Q for V = A (e 4,p + e QAlf ) (solid) vs. ACDM (dotted). The light 
yellow rectangles are the bounds < 0.1 from LSS, CMB, and BBN. 

The current investigation furthermore presents detailed simulations and 
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analyses of the evolution of the scalar field and its equation of state, the 
fractional energy densities in the scalar field, radiation, and matter, and the 
acceleration parameter in braneworld quintessential inflation for the sum of 
exponentials and cosh potentials. 



2 Cosmological Equations 

The homogeneous scalar field — since it is confined to the brane — still obeys 
the Klein-Gordon equation 

i + 3B*—% = -V.. (2) 

The Hubble parameter H is related to the scale factor a and the energy 
densities in matter, radiation, and the scalar field through the braneworld 
modified [4] Friedmann equation 



where the energy density 



P = P4> + Pm + Pr, P4> = \$ + V {(/)), (4) 

a is the four-dimensional brane tension, A 4 is the four-dimensional cosmolog- 
ical constant, and £ is a constant embodying the effects of bulk gravitons on 
the brane. We will set the four- dimensional cosmological constant to zero. 
The "dark radiation" term S/a A can be ignored here since it will rapidly 
go to zero during inflation. Thus for our purposes the modified Friedmann 
equation takes the form 

P (-x , P 



H = — ( 1 H — — I . (5) 

3MJ, \ 2a) V ' 

In the low-energy limit /) < <t, the Friedmann equation reduces to its stan- 
dard form H 2 = p/{ZM 2 P ). 

The conservation of energy equation for matter, radiation, and the scalar 
field is 

p + 3H(p + P) = (6) 
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where P is the pressure. Except near particle-antiparticle thresholds, P m = 
and P r = Pr/3. Equation (jBJ) gives the evolution of p m and p r , and the 
Klein-Gordon equation (j2J) for the weakly coupled scalar field, with 

P^ = \4> 2 -V{<p). (7) 

The acceleration equation 

follows from Eqs. (JSJ) and ([6]). While inflation occurs in the low-energy (stan- 
dard cosmology) limit when w = P/p < —1/3, for inflation to occur in the 
high-energy limit w < —2/3. 

We will use the logarithmic time variable r = ln(a/ao) = — ln(l+;z). Note 
that for de Sitter space r = H^t, where H\ = p\/(3Mp), and that H\t is a 
natural time variable for the era of A-matter domination (see e.g. Ref. |10j). 
For < z < zbbn ~ 10 10 , -23.03 < r < 0, and for < z < z Pl ~ 2.8 x 10 31 , 
-72.41 < r < 0. 

In Eq. (j3J), we will make the simple approximations 

p r = p r0 e" 4T , p m = p m0 e" 3r . (9) 

For the spatially homogeneous scalar field <f), the equation of state pa- 
rameter Wj, = w<f,(z) = P^/p^. The recent average of is defined as 

1 f T 

W = - w^dr. (10) 
r Jo 

We will take the upper limit of integration r to correspond to z = 1.75. The 
SNe la observations [11] bound the recent average w Q < —0.76 (95% CL) 
assuming Wq > —1, and measure the transition redshift z t = 0.46 ±0.13 from 
deceleration to acceleration (it is probably premature at this point to say 
more than that z t ~ 1). 

For numerical simulations, the cosmological equations should be put into 
a scaled, dimensionless form. Equations ([2]) and ^ can be cast [T5] in the 
form of a system of two first-order equations in r plus a scaled version of H : 

Hip' = ip (11) 
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(12) 




(13) 



(14) 



wh ere V = e 2T <p/H , H = e 2r H / Hq, V = e^V/p c0 , % = e 4T V v /p c0 , p = 
e 4r p/p c o, Pm = e* T pm/pco = ^moe r , pr = e 4T p r /p c0 = ^ro, <? = e 47 ~<* j ' p c0 , and 
where a prime denotes differentiation with respect to r: ip' = dip/dr, etc. 

This scaling results in a set of equations that is numerically more robust, 
especially before the time of BBN. 



Figure 2: Logio of H (blue, bottom) and H/H (cyan, top) vs. log 10 (l + z) 



Figure [2] illustrates that while H spans only ten orders of magnitude 
between ~ 10 23 and the present, H/H spans more than fifty orders of 
magnitude. 

3 Evolution of the Braneworld Universe 

In this section, we will closely follow the analyses of Refs. [5] and [I] . 
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3.1 Slow-Roll Braneworld Inflation 

The inflationary slow-roll parameter e is given by [I] 

H M 2 P fV^\ 2 1+V/a 



H 2 2 \VJ (l + V/(2a)) 2 
(The slow-roll parameter 



(15) 



^fiTW(*r £ (16) 

for the potentials considered here.) Inflation occurs for e < 1. The slow-roll 
parameter can be approximated during inflation by 

e « — (17) 

for the sum of exponentials potential, the simple exponential potential, and 
the cosh potential, since V ^> a and if ^> 1 for our models during inflation. 
Inflation ends when e = 1, implying [5] that the potential V e at the end of 
inflation is 

V e » 2AV. (18) 

To evaluate 14 and a, the COBE-measured amplitude of primordial density 
perturbations is matched against the theoretical value at N = 50 e-folds from 
the end of inflation, where 



yielding Vn ~ (iV + 1)14. Now matching against the amplitude [4J of density 
fluctuations 

(2xKT 5 ) = A|«^^^(l + ^) « go ^ 2 M ^ AV3 (20) 
gives [5] 

I*; 1 /* ~ i.n x io 15 GeV/A, a 1/4 « 9.37 x 10 14 GeV/A 3/2 . (21) 
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3.2 Gravitational Reheating and Initial Conditions 



At the end of inflation, gravitational particle production [13] [15] results in 
a radiation energy density [5] 

^ « 0.01^ « 2.8 x « 4.86 x lO" 17 ^ (22) 

where g p ~ 10-100 is the number of particle species that are gravitationally 
produced. This relation fixes the "initial conditions" redshift Z{ ~ 3.89 x 
lO 23 ^/ 4 / A at which gravitational particle production occurs, from the scaling 

Pr « 4.86 x lO- 17 (? p p « g p ^—9^L\ ~ (i + ^) 4 Pr0 (23) 
corresponding [5j to a reheating temperature 

where g>* (T) is the effective number of massless degrees of freedom at temper- 
ature T. The universe undergoes [5] a transition from the era dominated by 
the scalar field potential energy to an era of kination (p^ ~ 1/a 6 ) dominated 
by the scalar field kinetic energy as the p 2 /(2a) term in H 2 becomes small 
compared with p. For this era of kination to occur, the potential V must be 
sufficiently steep. Since during kination p T j ' p^ ~ a 2 , the universe eventually 
makes a transition [5] to the standard radiation dominated era, around a 
temperature T w 10 3 GeV/A, well before big-bang nucleosynthesis at T w 1 
MeV. For the sum of exponentials or cosh potential, the universe then evolves 
according to the QCDM model. The transition to an accelerating universe 
due to quintessence occurs late in the matter dominated era near z ~ 1, as 
in ACDM. 

The initial conditions at the end of inflation will be set in the high-energy, 
slow-roll limit. The initial value for fa is specified by Eq. (fl8|) . The initial 
value for fa follows from the high-energy, slow-roll limits 

n 2 V 2 

H 2 « — V- » -V- (25) 

and 

^^M^-^-^ < 26) 
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1 12 Ve 1,2 
-0 2 W — « -AV. 

2^ 6 3 



(27) 



4 Simulations 

For the computations below, we will use Eqs. (TTTj) (TT4T) with initial conditions 
specified at 2$ by and ^ oc ^ (Eqs. ([TBI and fl27j) ). We set A = 5 and g p 
= 100. The constant A in the potentials is adjusted so that Q^q = 0.74. This 
involves the usual single fine tuning. The final time tf is set by Zf = —0.999 
corresponding to tf ~ 100 Gyr. 

First we briefly summarize the properties of the exponential potential, 
and then turn our attention to the quintessential inflation potentials. 



V(ip)/A 


A/pco 


Pi 




w 


2 cosh(5<^) 


0.4 


48.0 


-0.03 


-0.87 


g5v _|_ gV^V 


5.7 


47.4 


-1.60 


-0.78 


e 5 ^ + 


2.2 


47.6 


-1.18 


-0.89 



Table 1: Simulation results for the cosh and sum of exponentials potentials. 



4.1 Exponential Potential 

The exponential potential V(ip) = Ae Xlp [T5] [T8] can be derived from M- 
theory [19J or from N = 2, 4D gauged supergravity [20J. 

For the exponential potential with A 2 > 3, the cosmological equations 
have a global attractor with = 3/A 2 during the matter dominated era 
(during which w<f, — 0); and with A 2 > 4, the cosmological equations have a 
global attractor with = 4/ A 2 during the radiation dominated era (during 
which u>0 = 1/3). For A 2 < 3, the cosmological equations have a late time 
attractor with = 1 and = A 2 /3 — 1. 

For A = \/2 and p m = 0, a — > asymptotically; if p m > 0, the universe 
eventually enters a future epoch of deceleration. In either case, there is no 
event horizon. For A < v^2, the universe enters a period of eternal accelera- 
tion with an event horizon. For A > \/2, the universe eventually decelerates 
and there is no event horizon. 
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The ACDM cosmology is approached for A < 1/ y/3. Significant acceler- 
ation occurs only for A ~ y/3. For A = a/3, Wo is much too high; for a viable 
present-day QCDM cosmology, A < y/2 [12] in the exponential potential. 

4.2 Cosh Potential 

In the cosh potential model 

V = 2 A cosh(5^) « p A cosh(5</?), (28) 

dark energy derives from the value of the potential near its minimum. This 
is the simplest way to "correct" the exponential potential to incorporate 
quintessence. 

In the simulations for the cosh potential (Figs. EHZJ), initially (p ~ 50 and 
then evolves to <po = —0.03 at to- The linear decay of <p(r) in Fig. 0] during 
the kination era occurs because p$ ~ |0 2 ~ 1/a 6 and thus <p' is a (negative) 
constant since both H and ip are proportional to e~ T in Eq. ( TlTi) . 

1 

0.8 
0.6 

G 

0.4 
0.2 
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log ( 1+z ) 

Figure 3: O for V = 2Acosh(5<^) (solid) vs. ACDM (dotted). 

Between z ~ 2.5 x 10 14 (as the universe enters its radiation dominated 
stage) and z « 10 (when « A 2 V/M|, becomes of order if 2 « p m /(3Mf,)), 
"sits and waits" at a small negative value. 
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Figure 5: Oscillating behavior of (p near z = for V = 2A cosh(5<£>). 
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Figure 6: w$ for V = 2Acosh(5<^). 




With the cosh potential, (p is very slowly oscillating and decaying about 
if = at the minimum of the potential for z ~ 1 (we have neglected any 
present-day particle production by the oscillating scalar field), with angular 
frequency d « \f?>\Hx and decay time constant t ~ 2Hj l 1 /3 (see Fig. [5]). 

Note that — > — 1 for t > £ (Fig- The acceleration parameter — g 
in Fig. [7] closely mimics the ACDM curve. 

4.3 Sum of Exponentials Potential 

In the sum of exponentials potential JT]), a shift ip — > </j + x simply results in 
a redefinition of the constants A and B: 

V{<p) ^V(<p + X )=A (e x<p + Be^) , .4 = e Xx A, B = e^ X)x B. (29) 

In a realistic particle theory A and I? would be set from first principles; here 
we choose B — 1. It will turn out then that A ~ p c o for quintessential 
inflation. 

To produce primordial inflation in the braneworld scenario with the sum 
of exponentials potential and the subsequent transition to the standard radia- 
tion dominated universe satisfying the constraints on Q^, X ~ 4.7 (this agrees 
with the estimate A ~ 4.5 given in Ref. [S]); for present-day quintessence, 
fi < \/2 [12]. We will choose A = 5 and /i = y/2 or 1 (these values appear 
naturally in low-energy limits of M/string theory): 

V((p) = A (e^ + e^) (30) 

or 

V(<p) = A (e 5 ^ + e*) . (31) 

For the sum of exponentials potentials (Figs. I8 HT51) . initially (p w 50 and 
then evolves to (po ~ —1. Again note the the linear decay of v 9 ( r ) i n Figs. [9] 
and [T3] during the kination era. 

Between z ~ 2.5 x 10 14 and z ~ 2 (when now m\ ~ fi 2 V/Mp becomes 
of order H 2 w p m / (3Mp)), (p "sits and waits" at a small negative value, and 
then becomes more negative. The linear decay of <f>{r) after to occurs because 
(p is approaching the late time attractor with = 1 and w = /i 2 /3 — 1 
(Figs. [TU1 and [T4|) . Whenever w = const and H approximates its standard 
form, ip(r) is linear in r since both H and if) are proportional to e ( - 1 ~ 3w ^ T / 2 in 
Eq. ffTTj) (ty = 1 during kination). 
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Figure 8: for 1/ = A (e 5 ^ + e^) (solid) vs. ACDM (dotted). 
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Figure 9: (f for V = A (e 5 ^ + e^) . 
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Figure 11: Acceleration parameter — q for V = A (e 5v + (solid) vs. 

ACDM (dotted). 
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Figure 12: Q for V — A (e 5 ^ + e v ) (solid) vs. ACDM (dotted). 
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Figure 13: tp for 1/ = A (e 5 ^ + e v ). 
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Figure 14: hi V = A (e 5 ^ + e*). 
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Figure 15: Acceleration parameter — q for V = A (e 5ip + e^) (solid) vs. ACDM 
(dotted). 



17 



Note the differences in the acceleration parameter — q in Figs. [TT] and [TS] 
for the two cases \i = \/2 (eventual deceleration with no event horizon) and 
H = 1 (eternal acceleration with an event horizon). 

5 Conclusion 

In the braneworld framework, the sum of exponentials and cosh potentials 
yield natural quintessential inflation scenarios, motivated by string/M the- 
ory. The quintessential braneworld models also provide — through an era of 
kination — a natural mechanism for the transition of the universe from the 
primordial ^-dominated era to the era of radiation domination. 

From the time of radiation domination to the distant future, mimics 

for the quintessential inflation potentials presented here. In the ACDM 
model and in braneworld quintessential inflation from a cosh or sum of ex- 
ponentials potential, there is a period between roughly 3.5 Gyr and 20 Gyr 
after the big bang when 0.1 < < 0.9. 

Note that for the cosh and sum of exponentials potentials Wq < —0.78 
(see Table 1), satisfying the observational bound w < —0.76. 

Braneworld quintessential inflation can occur for potentials with (cosh) 
or without (sum of exponentials) a minimum, and with (sum of exponentials 
with ii — 1 or cosh) or without (sum of exponentials with fi = \/2) eternal 
acceleration and an event horizon. 

In both the cosh and sum of exponentials potentials considered here, 
the low z behavior of provides a clear observable signal distinguishing 
quintessence from a cosmological constant. 

References 

[1] P. Binetruy, C. Deffayet , and D. Langlois, Nucl. Phys. B 565, 269 (2000) 
[arXiv:hep-th/99050T2l . 

[2] P. Binetruy, C. Deffayet, U. Ellwanger, and D. Langlois, Phys. Lett. B 
477, 285 (2000) [arXiv:hep-th/9910219| . 

[3] T. Shir omizu, K. i. Maeda, and M. Sasaki, Phys. Rev. D 62, 024012 
(2000) [arXiv:gr-qc/9910076] . 



18 



[4] R. Maartens, D. Wands, B. A. Bassett, and I. Heard, Phys. Rev. D 62, 
041301 (2000) |arXiv:hep-ph/9912464| . 

[5] E. J. C opeland, A. R. Liddle, an d J. E. Lidsey, Phys. Rev. D 64, 023509 
(2001) l[arXiv:astro-ph/000642l] . 



[6] D. N. Spergel et al, |arXiv:astro-ph/060 3449. 



[7] T. Bar reiro, E. J. Copeland, and N. J. Nunes, Phys. Rev. D 61, 127301 
(2000) |arXiv:astro-ph/9910214] . 

[8] A. S. Majumdar, Phys. Rev. D 64, 083503 (2001) 
|arXiv: astro-ph/0105518| . 

[9] N. J. Nunes and E. J. Copeland, Phys. Rev. D 66, 043524 (2002) 
[arXiv: astro-ph/ 02 04715] . 

[10] C. L. Gardner, Phys. Rev. D 68, 043513 (2003) 
[arXiv: astro-ph/0305080| . 

[11] A. G. Riess et al. [Supernova Search Team Collaboration], Astrophys. 
J. 607, 665 (2004) [arXiv:astro-ph/0402512| . 

[12] C. L. Gardner, Nucl. Phys. B 707, 278 (2005) |arXiv:astro-ph/0407604] . 

[13] L. H. Ford, Phys. Rev. D 35, 2955 (1987). 



[14] B. Spokoiny, Phys. Lett. B 315, 40 (1993) |arXiv:gr- qc/9306008|. 

[15] P. G. Ferreira and M. Joyce, Phys. Rev. D 58, 023503 (1998) 
|arXiv:astro-ph/9711102] . 

[16] C. Wetterich, Nucl. Phys. B 302, 668 (1988). 

[17] E. J. Copeland, A. R. Liddle, and D. Wands, Phys. Rev. D 57, 4686 



(1998) |arXi v:gr-qc/971 1068|. 
[18] M. Doran and C. Wetterich, Nucl. Phys. Proc. Suppl. 124, 57 (2003) 



[arXiv:astro-ph7 0205267] . 
[19] P. K. Townsend, JHEP 0111, 042 (2001) |arXiv:hep-th/0110072| . 



19 



[20] L. Andrianopoli, M. Bertolini, A. Ceresole, R. D'Auria, S. Fer- 
rara, P. Fre, and T. Magri, J. Geom. Phys. 23, 111 (1997) 
[arXiv:hep-th/9605032] . 



20 



